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Abstract 

This paper concerns with existence, uniqueness and asymptotic behavior of the solutions for a nonlocal 
coupled system of reaction-diffusion. We prove the existence and uniqueness of weak solutions by the 
Faedo-Galerkin method and exponential decay of solutions by the classic energy method. We improve 
the results obtained by Chipot-Lovato and Menezes for coupled systems. A numerical scheme is 
presented. 
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1 Introduction 

We consider the reaction-difFusion coupled system in parallel way via parameter a > of the form 



ut-a{l{u))Au + f{u-v) = a{u~v) in Q x {0, T) (1.1) 

Vt-a{l{v))Av-f{u-v) = a{v-u) in n x {0, T) (1.2) 

u = v = in dnx{0,T) (1.3) 

u{x, 0) = uo{x) in (1.4) 

v{x, 0) = vo{x) in fl (1-5) 



where u — u{x, t) and v = v(x, t) are real valued functions, is a bounded domain of K", dVl is 
the boundary of of class C^. / : R ^ R and a : R ^ M are Lipschitz continuous functions with 
a(^) > m > 0. / : L^(n) ^ R is a continuous linear form. 

For the last several decades, various types of equations have been employed as some mathematical model 
describing physical, chemical, biological and ecological systems. Among them, the most successful and 
crucial one is the following model of scmilinear parabolic partial differential equation, called the reaction- 
diffusion system 

Fill 

^-DAu^f{u) = 0, (1.6) 

where / : R" — > R" is a nonlinear function, and D is an n x n real matrix of diffusion. This reaction- 
diffusion model is obtained by combining the system of nonlinear ordinary differential equations called 
the reaction system 
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and the system of linear partial differential equation called the diffusion system 

^-DAu^O. (1.8) 

In 1998, L. A. F. Oliveira [T^ considered the reaction-diffusion system where D was a. n x n real matrix 
and / : R" M" is a function. He studied the exponential decay for some cases. Except for some 
pubhcations on the subject, such as the searching for travehng waves solutions and some problem in 
ecology and epidemic theory, most of authors assume that diffusion matrix D is diagonal, so that the 
coupling between the equations are present only on the nonlinearity of the reaction term /. However, 
cross-diffusion phenomena are not uncommon (see [1] and references therein) can be treated as equations 
like in which D is not even diagonalizable. In 1997, M. Chipot and B. Lovat [5] studied the existence 
and uniqueness of the solutions for non local problems 



ut-a\^J u{x,t)dxj Au ^ f{x,t) in f7 x (0, T) (1.9) 

u{x, t) ^ on dVL X (0, T) (1.10) 
u{x, 0) = uq{x) on Q. (1-11) 

where £7 is a bounded open subset in K" n > 1 with smooth boundary d^. T is some arbitrary time, a 
is some function from R into (0, -l-oo). This problem arises in various situations, for instance u could de- 
scribe the density of a population(for instance of bacteria) subject to spreading. The diffusion coefficient 
a is then supposed to depend on the entire population in the domain rather than on the local density i. 
e. moves are guided by considering the global state of the medium. They proven the following result: 



Theorem 1.1. Let Tq > 0, uq G ^^(^)i '^o > 0, Uq ^ 0. Let a be a continuous function positive in 
a neighborhood of u^dx. Then for f G i^([0, T] : H^^{VL)) there exists < T < Tq and u solution to 
(f01l - (frTT|) such that 



ueL\[0, T] : i/i(f7))nC"([0, T] : L^n)), 
GL'([0, T] : H-\n)), 

<ut, V > + a(^J udx^ {\/u-Vv) ^< f, V > V w G iJp (f7), a.e. t e [0, Tq] 
where (Vu • Vw) = Vm • Vu dx. 

In 2005, S. D. Menezes ff, give a simple extension of the result obtained by M. Chipot and B. Lo- 
vat [5], considering a = a{l{u)), f = f{u) continuous functions. Indeed, they studied the existence, 
uniqueness and periodic solution for the following parabolic problem 

ut -a{l{u)) An + f{u) = /i in f7 x (0, T) (1.12) 
u{x, t) = on dn X (0, T) (1.13) 
u{x, 0) = uo{x) on Q (1-14) 

where is a bounded open subset in R", n > 1 with smooth boundary dfl. T is some arbitrary time. 
I : L^{^) ^ R is a nonlinear form, /i e L^(0, T : H^^{n)) and T > is some fixed time. This problem is 
nonlocal in the sense that the diffusion coefhcient is determined by a global quantity. This kind of prob- 
lems, besides its mathematical motivation because of presence of the nonlocal term a{l{u)), arises from 
physical situations related to migration of a population of bacteria in a container in which the velocity of 
migration 1? = a Vw depends on the global population in a subdomain fi' C fi given by a = ^(/q/ udx). 
For more information [S] and reference therein. 

This article is concerned with to prove the existence, uniqueness and the exponential decay of the system 
(ll.ip - (|1.5p using the energy method. The method of energy consists of to use appropriate multipliers 
to build a functional of Lyapunov, in this direction we prove that for this types of materials where the 
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energy, that can flow from one part to another, is strong enough to produce exponential decay for the 
solution of the coupled system. 

This paper is organized as follows. Before the main result, in section 2 we briefly outline the nota- 
tion and terminology to be used subsequently. In the section three we prove the existence and uniqueness 
of solution, in the section four we prove the exponential decay of solution of the system. Finally, numer- 
ical evidence corroborating our theoretical results is given in section five. In this paper, we prove the 
following two theorem: 

Theorem 1.2. Let (uq, vq) €E L^(ri) x L^(0) and < T < +oo, where the time T depends only |mo|l2(s^) 
and |?;o|l2(o). // (fTHll- (fT^ holds, then there is at most one solution of (frT|) - (fT3|l in L'^{0, T : H^{n))n 
C([0, T) : L^in)) X L2(0, T : H^{n)) n C{[0, T) : L^{n)) with initial data {u{x, 0), v{x, 0)) = {uq, vq). 

Theorem 1.3. Let (u, v) be a solution of system (ll.ip - dl.Sp given by the theorem L2, then there exist 
positives constants C and rj, such that 

E{t) <CE{0)e-"*. (1.15) 



2 Preliminaries 

In this work we consider the reaction-diffusion coupled system in parallel way via parameter a > as 

ut~a{l{u))Au + f{u-v) = a{u-v) in Q x {0, T) (2.1) 

Vt-a{l{v))Av-f{u-v) = a{v-u) in Q x {0, T) (2.2) 

u = v = in dflx (0, T) (2.3) 

u{x, 0) = uo{x) in n (2.4) 

v{x, 0) = vo{x) in n (2.5) 

where is a bounded domain of M", dil is the boundary of fl of class C^. / : R ^ R is a Lipschitz 
continuous function, that is, there exists Mi > such that 

\fis)- f{t)\<Mi\s-tl Vs, teK. (2.6) 

a : R ^ M is a Lipschitz continuous function, that is, there exists A/2 > such that 

\a{s)-a{t)\<M2\s-t\, V s, t G K. (2.7) 

with 

a(O>TO>0, VCeM (2.8) 

and 

I : L'^{Q) M. is a continuous linear form. (2-9) 

In the system, the distributed spring coefficient is coupled by the terms a (u — v) and a {v — u). In this 
sense the Energy can flow from one part to another through this parameter a. 

By < ■ , • > we will represent the duality pairing between X and X' , X' being the topological dual 
of the space X, We represent by H"^{^1) the usual Sobolev space of order m, by i/™(il) the closure of 
C^{Vl) in i/™(il), and by L^(ri) the class of square Lebesgue integrable real functions. In particular, 
(r2) has inner product (( • , • )) and norm || • || given by 



{{u, v)) = / Vu-^vdx and = / |Vupda; 

Jo. Jn 
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For the Hilbert space L^(rj) we represent its inner and norm, respectively, by ( • , • ) and | • |, defined by 

{u, v) — / uv dx and |wp = I dx. 

n 



Throughout this paper c is a generic constant, not necessarily the same at each occasion(it will change 
from line to line), which depends in an increasing way on the indicated quantities. 

We take the initial conditions as following 

{uo{x), vo{x)) £ L^{n) X L^{n). (2.10) 
We denote the potential energy associated to this system by 

m = ljj\u\' + \v\'] dx. (2.11) 



3 Existence and Uniqueness of a local solution 

In this section, we will prove that for {uq, vq) G L^(ri) x i^(ri) there exists a unique solution of (|2.6p - 
(HH) in L2(0, T : H^in)) n C([0, T) : L^in)) X L^{0, T : H^ifl)) n C([0, T) : L^in)) where the time 
T depends only iMoli^^fi) and |wo|L2(n). We make use of Faedo-Galerkin approximation for to prove the 
existence of weakly solutions. We write the system p.l|) - (|2.5ll in the following form. 



Ut-a{l{u))Au = - f{u, v) + g{u, v) in Q = n x {0, T) (3.1) 

vt-a{l{v))Av = f{u, v) ~ g{u, v) in Q = n x {0, T) (3.2) 

u = v = on dnx{0,T) (3.3) 

u{x, 0) = uo{x) in n (3.4) 

v{x, 0) = vq{x) in fl (3.5) 

where we denote f{u, v) = fiu — v) and g{u, v) = a{u ~ v). 

Theorem 3.1 (Existence). Let [uq, vq) e L^{n) x L^{n) and < T < +oo. // (lMl)-(l231) holds, then 
there exists (u, v) solution of (|3.ip - (|3.5p such that 

{u, v) G L^(0, T : H^(n)) n C([0, T) : L^{n)) x ^^(o, T : H^{n)) D C{[0, T) : L^{n)) (3.6) 

{uu vt) e L^{0, T : H-\n)) x L^iO, T : H-^Q)) (3.7) 

— {u, hi) + a{l{u)) [ Wu-Whidx — — f f{u,v)hidx+ f g{u,v)hidx (3.8) 
dt Jn Jn Jn 

for all hi e i?g(ri), where p.Sp must be understood as an equality in X''(0, T). 

— {v,h2)+a{l(v)) j Vw • V/i2 = / f(u,v)h2dx~ f g{u,v)h2dx (3.9) 



dt 

for all /i2 G HQ{fl), where p.Op must be understood as an equality in T>'{0, T). 

Proof, i) Approximate problem: Let {uijjjgN be a Hilbertian basis of Hl{Vl) (cf. H. Brezis, [3]). 
For each j = 1, 2, 3, . . . represent by Vj, the subspace of i?o(f^) generated by {wi , W2, • ■ • , Wj}. The 
approximate problem, associated with (|3.ip - l|3.5p . consists of to find itj, Vj S Vj such that 

{u'j, hi) - a{l{uj)){Auj, hi) = - {f{uj, vj), hi) + {g{uj, Vj), hi), y hi e Wj (3.10) 

{v'j, h2) - a{l{vj)){Avj, /12) = if{uj, Vj), ^,2) - {g{uj, vj), ^-2), V /12 e Vj (3.11) 

Uj{0) — uoj Uq strongly in L'^(f2) (3.12) 

Vj{0) = voj Vq strongly in i^(57) (3.13) 
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Let hi = Wi{x) and /12 — Wi{x) for i = 1, . . . , j, then in (|3.10|) - (|3.13|) we have for 9kj, <j>kj £ C°°{^1) 



\fc=i / \ \fe=i / 

+ [s [ '^dkj{t)wk{x), '^(f>kj{t)Wk{x) , Wi{x) 



^'^dkj{t)wk{x), Wi{x) 



k=l 



k=l 



and 



^<l^'kji't)wk{x), w,{x) 1 - « I M ^(t)kj{t)wkia 

k=l / \ \k=l 

f ^9kjit)'Wkix), ^(l)kjit)wkix) , Wi{x) 



A'^(f)kj{t)wk{x), Wi{x) 



k=l 



\k=l 



k=l 



9 '^dkj{t)wk{x), '^(l)kj{t)wk{x) , Wi{x) 
\ \k=i k=i J / 



that is, 



^kji*) - >^ka{l{uj))ekj{t) = - {f{u.j,v.j),Wi) + {g{uj,vj),Wi) 
^'kj{t) - Afe a{l{vj)) (t>kj{t) = ifiuj, Vj), w,) - {g{uj, Vj), w,). 



(3.14) 
(3.15) 



ii) Approximate solutions: We will just work with the equation p.l4p . For (|3.15p the result is similar. 
For i , A; = 1 , . . . , j in p.l4p , we have the following system 



9' 



JJ 



Ai a{l{uj)) 




that is, 



\j a{l{uj)) 



X' 
X{0) 







' (/, wi) - 




(.9, wi) 






(/, ^2) 


+ 


{9, W2) 










. (9, Wj) _ 



Fix, t) 



(3.16) 
(3.17) 



where F{X, t) = AX + B and Xq = [ay, a2j, ctjjY ■ The system p.l6p - p.l7p is equivalent to 
system of ordinary differential equations of first order. Let us show that the system p.l6p - p.l7p is in 
the conditions of Caratheodory's theorem. 

Claim. For fixed X, we will show that A and B are measurable in t. 

In fact, we observed that the matrix A is formed for the elements Xka{l{uj)) with k = 1, 2, . . . , j. 
Since / is a lineal and continuous form and the operator a is continuous, then the composition a{l{uj)) 
is also continuous; therefore a{l(uj)) is continuous for k — 1, 2, . . . , j and then A is measurable in t. 
On the other hand, let us observe that B is formed by the elements (/(itj, vj), Wi) and {g{uj, Vj), Wi), 
with i = 1, 2, . . . , j. Since / and g are continuous and Wi G ifg(ri), we concludes that B is continuous 
and therefore measurable. 



Claim. For fixed t, we will show that F is continuous in X. 
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In fact, notice that B is continuous in X, because B is constant in relation to X. For continuity of 
AX, is enough we verify that A is continuous in X. Let Yik — ^kj {k — 1, 2, . . . , j) be the projection 
W — > M and a{X) = Yl^ X Wk- For each t fixed, as 

j 

'^jW = y^Jkjjt) Wk, 

k=l 

we can consider the function 

X ^a{liu,))=a(^l (j29k,it)wk)i^ = « 11 ^ ) ' 

Since A is Uneal combination of continuous functions, proceeds that A is continuous in X, hence the 
function F(X, t) is continuous in X. 

Claim. Let if be a compact of D = E x [0, T], where E ^ {X G : ||X||r,xi < S, 6 > 0}. 

We will show that exists a real function mr{t), integrable in [0, T], so that 

\\F{X,t)\\R,.i <mr{t), y{X,t)eD. 

In fact, we denote by || • \\pq the norm of maximum in M.p'^. But F{X, t) = AX + B, then 

m,xi < iixii,xi + iii?ii,xi. 

Since X £ E, we have \ \X\\jy,i<5. Then 

||^^(^, t)||,xi <5||A||,x, + ||i?||,xi. 
Notice that Xk a{l{uj)) are continuous functions, hence ||^||jxj < C (C > 0). 
On the other hand, for the matrix B we have 

\{fiUj, Vj), W^)\ < \f{Uj, Vj)\ \w^\ = \f{Uj, Vj)\, 

\igiuj, Vj), w^)\ < \g{uj, Vj)\ \w,\ = \g{uj, Vj)\. 
Thus i)||jxi < S C + \f{uj, Vj)\ + \g{uj, Vj)\ = mr{t), where mr{t) is integrable in [0, T]. 

Hence, the system p.l6p - p.l7l) satisfies the conditions of Caratheodory, and then exists {uj{t), Vj{t)} e 
[0, tj) X [0, i,), t, < To. 

We now have to establish an estimate that permits to extend the solution {uj{t), Vj{t)} to the whole 
interval [0, T]. 

From now on, {Ci}i=i...7, will denote positive constants, independents of j and t. 

iii) A priori estimates: We put hi — Uj and h2 = Vj in the equations (|3.10p and (j3.1ip respectively, 
we have 

{u'j, Uj) - a{l{uj)) {Auj, Uj) = - {f{uj, Vj), uj) + {g{uj, Vj), Uj) (3.18) 
(Vj, vj) - a{l{vj)) {Avj, Vj) = {f{uj, Vj), Uj) - {g{uj, vj), uj). (3.19) 

Using the boundary condition and the first Green's identity we have 

{—Auj, Uj) = / {—Auj)uj dx — I \7uj ■ Vuj dx = iVwjp — 
Jfl Jfl 
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Then, we can write p.lSp as 
1 d 



2dt 



\uj\ + a{l{uj)) \\uj\\ ^- {f{uj,Vj),Uj) + {g{uj,Vj),Uj). 



In a similar way we can write p.l9p as 
1 d 



2 dt 



\vj\ +a{l{vj))\\vj\\ ^ {f{uj, Vj), Uj) - {g{uj,Vj), Uj). 



(3.20) 



(3.21) 



Adding the equations (|3.20p with (|3.2ip and integrating of to t, we obtain 



[a{liu,))\\u,\\^ + a{liv,))\\v,\\^]ds = 0, 



WM^ + + 2 / [a{l{u,)) \\u,\\^ + ail{v,))\\v,\\' 1 ds = \u,m + |t;,(0)|. 



that is, 



In the last identity, using (|2.6p and (|2.8p we obtain 

I^WP + hit)? + 2m [\ ll^.lp + \\v,\\^ ) ds < \u,m^ + \vM?- 
Jo 

Since Uj{0) uq and Vj{0) vq strongly in i^(il) follows that |uj(0)|^ + |wj(0)|^ < C. Hence 

\u,it)\' + \vM' <c. 
From where follows that Uj{t) and Vj{t) are bounded in L°°{0, T : _L^(il)). Thus, 

(iKip + ii«,inds<c, 

then Uj{t) and Vj{t) are hmited in L^{0, T : H^{Q)). 

From ([XTU1) - ([XTT|1 . we have that 

u'j — a{l{uj))Auj — /(wj, Vj) + g{uj, vj) e H^^{fl) 
= a{l{vj))Avj + f{uj, Vj) - g{uj, v^) £ H-^{n). 

Notice that — a{l{uj)) Auj defines an element of H^^{^1), given by the duality 

{-a{l{uj))Auj, hi) ^a{l{uj)) Vuj-Vhidx, V /ii £ iJg (f^). 

Jn 

In a similar way we have 

{-a{l{vj))Avj,h2)^a{l{vj)) Vuj • V/i2 dx, V /i2 S i?o (f^). 



(3.22) 



Using the fact that — a{l{uj)) Auj, — a{l{vj)) Avj G H ^(fi), the dual of iJg (fi), then they are lineal and 
continuous forms and therefore both are bounded. 

Since uj, Vj e L^{0, T : L'^ifl)), then 



\f{uj,Vj)\dx< / (3\uj — Vj\dx< / (|uj| + jujl) 
In Jn 

1/2 / » ^ 1/2 

'n 



< C 



In '"^ '^'^^) 



\vj\^dx 



and 



\g{uj , Vj)\ dx — / a\uj~Vj\dx< / a {\uj\ + \vj\) dx 

Jn 

1/2 / /. \ 1/2" 



< c 

Therefore 



J \uj\^ dx^ + (^j \vj\^ dx 



and we concludes that u^-, v'^ are bounded in L^{0, T : H^^{Vl)). 

iv) Passage to the limit: We have that 

Uj, vj are bounded in L°°(0, T : L^ifl)) n ^^(O, T : H^ifl)), (3.23) 



u'-, w:,- are bounded in L^{0, T : H-\n)). (3.24) 



4' 



Now, due the coroUary of Banach-Afouglu (See [13], p. 66), we can extract subsequences of uj^ = uj and 
Vj^ = Vj (which we denote with the same symbol) so that 

Uj ^ u weak star in X°°(0, T : L^{fl)) (3.25) 

Vj V weak star in i°°(0, T : L^{fl)) (3.26) 

Uj ^ u weak in L'^{0, T : H^{fl)) (3.27) 

Vj V weak in L^{0, T : H^iQ)). (3.28) 



Consequently 



{uj, hi)dt — > / (m, hi)dt, V/ii e L°°(0, T : L^{n)) (3.29) 
Jo 

(uj, hi) dt — > [ {u, hi) dt, VhiE L^{0, T : H^{n)) (3.30) 
Jo 

T rT 



ivj,h2)dt — > {v,h2)dt, y h2e L^iO,T : L'^in)) (3.31) 
Jo 

[ {vj,h2)dt — >[ {v,h2)dt, V h2e L\0,T : H^{n)) (3.32) 
Jn Jo 



For (jX^ it proceeds 



^ u weakly in ^^(O, T : H-^{n)) (3.33) 



4 ^- 



weakly in L^{0, T : H-\Q)). (3.34) 



On the other hand, i?o(f^) ^ i^(ri) ^ if ^i^)- By Lions- Aubin's compactness Theorem [TT] follows 
that 



— > u strongly in L^(0, T : L\n)) (3.35) 

— >v strongly in L^{0, T : L^{n)). (3.36) 

The convergence p.25l) - (l3.26l) means that 

i-T pT 

{uj{t), w{t))dt — > {u{t), w{t))dt, Vw eL\0,T : L^{n)) 



T 



{vj{t), w{t))dt — y {v{t), w{t))dt, yweL\0,T:L^{n)). 
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We choose w = 9hi with 6 £ V{0, T), hi S L'^{n) and we wiU show that for aU 9 e ^(O, T) and for aU 

hi e L^n), 

T 

[{g{uj, V,), hi) - {g{u, v), hi)]e{t) dt 0. 



Let T be a positive number such that supp{9) C [0, T], then 

[{giuj, Vj), hi) - {g{u, v), hi)]9{t)dt = / {g{uj, Vj) - g{u, v), hi)9{t)dt. 
Jo 

Hence, by straightforward calculations 

{g{uj, Vj) - g{u, v), hi)9{t)dt < \g{uj, Vj) - g{u, v)\\hi\\9{t)\ dx dt 

Jo Jn 

= / / \a {uj — Vj) — a {u — v)\ \hi\ \9{t)\ dx dt 
Jo Jn 

n\a{uj -u) ~a{vj ~v)\ \hi\ \9{t)\dxdt 

< [ [ {a\uj - u\+ a\vj - v\)\hi\\9{t)\dxdt. 
Jo Jn 

Using L^(0, T : L'^{^)) L^{0, T : L'^{fl)) and the Cauchy-Schwartz inequality we obtain 

^dxj U \hi\^dxj dt 



{g{uj, Vj) - g{u, v), hi)9{t)dt < C I \uj - u\ 
Jo \Jn 

T / I- ^ 1/2 / » \ 1/2 



C [ ([ \vj-v\'^dx\ ([ \hi\'^dx] dt. 
Jo \Jn J \Jn J 



Applying the Cauchy-Schwartz inequality and considering the convergence (I3.35P we obtain 

1/2 / „ X 1/2 



C f (f \uj-u\^dx\ (f \hifdx\ dt 
Jo \Jn J \Jn J 



|2 J„ jj- I I / / I J. |2, 



< / / \uj -ul-'dxdt] / / \hi\'dxdt\ <e. 
\Jo Jn J \Jo Jn J 

In a similar way using the convergence (|3.36p we have 

C [ ( [ \vj-v\^dx] ' ( [ \hi\'^dx] ^ dt 
Jo \Jn J \Jn J 

< I I \vj-v\^dxdt\ \ ( I \hi\^dxdt\ <e. 



Jn I \Jo Jn J 



Therefore we have 

fT 



{g{uj, Vj) - g{u, v), hi) 9{t) dt < e. 







Performing similar calculations we can to prove that 







(/(uj, v.j) - f{u, w), hi)9(t)dt < e. 



We will show now, that for every 9 E V ([0, T]) and for every hi € _L^(r2) 

a(^K)) / / Vuj ■Vhi6{t)dt — > a{l{u)) [ [ Vu ■ Vhi 6{t) dt. 
Jo Jn Jo Jn 



(3.37) 



It is enough we prove that 

a{l{uj)) — >a{l{u)) in £^(0, T), V T > 0. (3.38) 
Since a is continuous, we will show that 

l{uj) — >l{u) strongly in L'^{0,T). (3.39) 

In fact, because 



T 

2 



\l{uj)-l{u))Ydt^ \l{uj ~u)Ydt <Cq \uj-u\^dt<e. 
Jo Jq 

This last one result, is consequence of the convergence p.35p . 

These convergence implies that we can take limits in the approximate problem (|3.1ip - (|3.15p . and then to 
verify the conditions (i), {ii), [Hi) and (iv) of the Theorem. 

Now, we will make verify of the initial data and we prove the uniqueness of solutions. Using the re- 
sult of regularity we have that 

w, V e C"(0, T : L^{n)) (3.40) 

In this form, makes sense we calculate u{0) e w(0). Let us consider 9 E C^(0, T : M), with 9{0) = 1 and 
9{T) = 0. Since the convergence p.29p we have 

T pT 

2/ 



{u'j,z)9dt — >/ {u',z)9dt, zEL^{n). (3.41) 
Jo 

Performing integration by parts in (|3.4ip we have 

-{uj{0),z)-[ [uj,z)9'dt — >-{u{0),z)-[ {u,z)9'dt. (3.42) 



Using the convergence p.29p in (|3.42p we have (uj(0), z) — > (u(0), z), for all z E Hl(n). But Uj{0) 
converges strong for wq in L^(ri), consequently weak in L'^{Vl). Therefore (uj(0), z) — > (uq, z), for all 
z E Hliyi). From uniqueness of the limit we have (u(0), z) — > (uo, z), for all z E H}^{Vl). Thus, u{Q) — uq. 
In a similar way we can prove that v{0) = vq- 

To finish this section we will show the uniqueness of solution. 

Theorem 3.2(Uniqueness). Let {uq, Vq) E L^{il) x L'^{n) and < T < +oo, where the time T de- 
pends only |uo|l2(o) a?^^ I''^o|l2(s7) . // (|2.6p - (|2.9p holds, then there is at most one solution of (|3.ip - (|3.5p 
in L2(0, T : H^iQ)) n C([0, T) : L^{n)) x L^{0, T : H^{n)) D C([0, T) : L^iQ)) with initial data 

{u{x, 0), V{X, 0)) = {uq, Vq). 

Proof. Assume that (ui, vi), (wz, ^2) in L^iO, T : H^in)) n C([0, T) : L^{VL)) x L^{Q, T : H^{n)) n 
C([0, T) : L'^in)) are two solutions of (EUl-dXH) with ut, vt in L^{0, T : R-^Q)) x L^{0, T : H^'^{n)), 
so all integrations below are justified and with the same initial data, in fact (ui — U2){x, 0) = and 
(U2 - U2)ix, 0) = 0. Then 

■^ui - a{l{ui)) Aui = -f{ui,vi)+g{ui,vi) 
dt 

^vi - a{l{vi)) Avi = f{ui,vi)-g{ui,vi) 
at 
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and 



^U2 - a{l{U2)) AU2 = -f{U2,V2)+g{U2,V2) 
^V2 - a{l{v2)) AV2 = f{U2, V2) - g{U2, V2). 



Then 



— {ui—U2,hi)+a{l{ui)) f Vui ■ Vhi dx — a{l{u2)) f Vm2 • V/ii da; 

Jn Jn 

= -{f{ui,vi)-f{u2,V2),hi) + {g{ui,vi)-g{u2,V2),hi), V/iie-ffo(0) 

■^{vi - V2, h2)+a{l{vi)) [ Vvi ■S/h2dx - a{l{v2)) [ Vv2-Vh2dx 
at Jn Jn 

= (/(wi, vi) - f{u2, V2), h2) - {g{ui, vi) - g{u2, V2), /12), V /i2 e Hq{Q.). 
Using that g{ui, vi) — g{u2, V2) = 01 {u\ — U2) — a{v\ — V2) and that f{u, v) = f{u — v) we have 

— {ui—U2,hi) + a{l{ui)) I Vui ■ Vhidx ~ a{l{u2)) [ Vu2-Vhidx 
dt Jq Jn 

= -{f{ui-vi)-f{u2-V2),hi)+a{ui-U2,hi)-a{vi-V2,hi), V/iiGi/"o(^) 

■j:{vi — V2, /i2) + a{l{vi)) [ Vvi ■ V/12 dx — a{l{v2)) [ Vt;2 • V/12 dx 
dt Jn Jn 

= {f{ui-Vi)-f{U2-V2),h2)-a{ui-U2,h2)+a{vi-V2,h2), V/l2GiIo(^)- 

On the other hand, let hi = ui — U2 and h2 = v\— V2 we obtain 

^\ui — U2\'^ + a{l{ui)) [ ^ui ■V{ui — U2) dx — a{l{u2)) [ Vu2 ■ y{ui — U2) dx 
dt Jn Jn 

2 



and 



Hence 



and 



= -{f{ui-Vi)- f{U2-V2),Ui-U2) + a\ui-U2\ - a {vi - V2, Ui - U2) 



^\vi - V2\^ + a{l{vi)) [ Vvi-V{vi-V2)dx~ ail{v2)) [ Vv2 ■ V {vi - V2) dx 
dt Jn Jn 

= {f{ui - Vi) - f{U2 - V2), Vi -V2)-a (Ui - U2, Vi -V2) + a \vi - V2f. 



i\ui - U2f + o(«(«i)) + ail{u2)) \\u2\f + [ail{ui)) - a{l{u2) ] ( Vui • Vu2dx 

dt Jn 

< \f{ui - Vi) - f{u2 - V2)\ \ui - U2I + a \ui - U2\^ +a\vi- V2\ \ui - M2I 



^^\v^-V2\'' + a{l{v{)) \\vif + a{l{v2)) \\v2\\'' + [a{l{v^)) - a{l{v2)] j^Vvi ■ Vv2dx 

< \f{ui - Vi) - f{U2 - V2)\ \V1 -V2\+a \ui - U2I \V1 -V2\ + a \vi - V2\'^. 
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Using (|2.6p - (|2.9p and the Young inequality we have 
d 



dt 



|'"i~U2| +m||Mi|| +m||u2|| 

< Mi\l{ui)-liu2)\ \\U1\\\\U2\\ + A'hliui ~ U2) ~ (Vi - V2)\ \UI-U2\ 



I |2 I |2 I 

a\ui~U2\ + — \vi-V2\ + ^ l"l 



-"2 



< Ml Ci \ui ~ M2I ||ui|| ||"2|| + Afs |mi - + M'i \ux ~ U2I \vx - V2\ 

^1 1 2 ^ I i2 

+ 2" pi ~ "21 + l^l - V2\ 



m ,,n M?C?,, , ,9 ,^ , ,9-^3, 



W2 



^3 I ,2 



-a |ui - U2 



< 



IKI 



^^2 (72 3 

^ IIM2IP |ui - M2p + - (M3 + a) |ui - W2p + (A/3 + a) - wal 
2to 2 2 



and 



Then 



and 



< M2\l(yi)-l(y2) I llwill IIW2II +M3|(wi -U2) - (wi - W2)| |wi - W2I 
+ - |ui - U2p + |fi - W2p + a \v\ - W2p 

< Af2C2 \VX~V2\ \\vx\\ IIW2II +M3 |wi - W2p +M3 |mi - W2I ll^l - W2I 

|2 



f2 r<2 



"2 



< 



m ,, ,,2 , M2^C2" ,, ,,2, ^, ,2 , ,^ , , ^, ,2 , Af3 , ^, ,2 



2 11-11^ + 2^ 



\\V2\\ \vi - -^2! + M3 - V2Y + — |wi - ^2!' 



— \ux-U2\ 



|2 " I 
W2I + — |Wl - "2 



..2 , M2^C2^ ,, ,,2, ,2 3 ,,,, . , ,2 1 



2m 



I|w2ir K'l - f'2r + 2 (^'^3 + a) iwi - V2r + {M3 + «) im - U2I 



-M2|' + y ||"i||'+m||?/2||' 



2m 



< 



\\U2\\^ \ui -U2\^ + - {Ah + a) \ui - U2p + - (Af3 + a) \vi - V2[' 



(3.43) 



^|i'i-f2|' + y ||i;i|P+m||v2||' 



< 



2m 



\v2\\'^ \vi - W2p + 7: (Af3 + a) \vi - V2\^ + - (Ms + a) \ui -U2p 



(3.44) 



Adding ([3^ with (|04|) we obtain 



^( - W2p + ifi - «2n + Y ikiir 

A/?C?,, ,,2, ,2 M^Ci 

< ^^||"2||'|W1-W2|' ' ' ' 



-m\\u2\\^ + - 

|W2|P \vi - W2 



||wi|P + m||wi 



2m 2m 

+ 2 {Ms + a) \ui - U2p + 2 (M3 + a) |ui - t;2p 



2m 



||u2||^ + 2(Af3 + a) 



\ui - U2I 



M^|h|P + 2(M3 + a) 
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We define 



2m 2m 
Thus, 

^ [ |m1 - + K'l - t'2p] < (pCi) |W1 - + ~ f'2p- 

Let 7^(^) = sup{.^(i), ^(t)}, then 7^ > and 

^ ( |U1 - U2P + l^-l - ^21') < n{t) ( |ui - + l^'l - V2|^ ) . (3.45) 

Integrating p.45p over t e [0, T] and using that mi(0) = ^2(0) and ui(0) = ^2(0), we obtain 

|mi - ""2!^ + l^i - ^2!^ < / Hit) [\ui - U2\^ + \vi ~ V2\^) dx. 
Jo 

Let p{t) — \ui — W2p + \vi ~ ^2!^, then 

p{t) < [ n{s)p{s)ds. (3.46) 



Jo 

Applying GronwaU's inequahty, we obtain 

p{t) < 0. 

Therefore p = 0, i. e., 

\ui - -^2^ + \vi - t;2p = 0. 
Using the regularity of the solutions, the uniqueness follows. 

4 Exponential stability 

In this section we show that the total energy (|2.1ip associated to system (|2.ip - (|2.5p decay exponentially 
to zero as t tends to infinity. In what follows we will prove our main result: 

Theorem 4.1. Let (u, v) be a solution of system ()l.ip - (|1.5p given by the theorem 3.1 and theorem 
3.2. We suppose that m > 2cp(Mi + a) > 0, where Cp corresponds to the constant of the Poincare 
inequality.. Then there exist positives constants C and r/, such that 

E{t) <CE{0)e-"*. (4.1) 

Proof. Multiplying equation (|2.ip by u and integrating over x £ $7 we have 

1 d 
2di 



\u\ dx + a{£{u)) / |Vw| dx — — f{u — v)udx + a j {u~v)udx. 
n Jn Jn 



Multiplying equation (|2.2p by v and integrating over x D, we have 

-— f \v\'^ dx + a{i{v)) f \Vv\'^ dx — ( f{u — v)vdx — a ( {u — v)vdx. 
2 dt Jq Jq 

Adding the expressions above and using (|2.8p we have 

1 d 



2dt 



f + Iwp] dx + m / [|Vup + |Vup] da; < / \f{u-v)\\u-v\dx + a \u~v\'^dx. 

n Jn Ja Jn 
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Using (|2.6p follows that 
d 



^^E{t)<-ml [ |Vwr + |Vwr] da; + (A/i + a) / \u - dx. 



Applying the Poincare inequality and the Young inequality, we obtain 

d 9m 

^E{t) < Eit) + 4 (Afi + a) E{t). 

dt Cp 

Now choosing m > 2cp {Mi + a) > follows that there exists C > such that 

±E{t)<-CEit). 

From we concludes that there exists 77 > such that 

E{t) < CE{0)e-'^K 

The proof follows. 

5 Numerical Results 

In this section we consider a particular case of the nonlocal reaction diffusion equations (|2.ip - (|2.5p in 
one dimensional space {n — 1), fl = (0,1), and / : L^{fl) R defined by l{u) — / u(x)dx. Then we 

approximate the solution of the system using implicit finite differences. The numerical scheme reads as 
following 

+ /K - <) = a[u\ - <), (5.1) 




-/(u^-«f) = a(«f-«f), (5.2) 

(5.3) 

for i = 1, . . . , J — 1 and A; = 0, . . . , if, where 8t ~ T / 5x = 1/ J ^ Xi ^ iSx, i — 0, . . . , J and tk — k6t, 
k = 0, . . . ,K. In ()5.ip - (|5.3p . u'^ denotes the approximation of u{tk,Xi). In order to solve the system 
()5.ip - (|5.3p . we consider the initial condition approximated by 

u^=UQ{xi), and v^—VQ{xi), for i — 0, . . . , J. (5.4) 

For each k — 0, . . . , K, the scheme ()5.ip - (|5.3p is equivalent to a linear system with a tridiagonal matrix 
of M^'^"-'^)^'^''"^^ which is positive definited and then there exists a unique solution of ()5.ip - ()5.4p . 

In order to compare the numerical behaviour of the solution with theoretical and numerical behaviour 
of the solution for one equation (scalar case), we take the same nonlinear reaction terms with similar 
parameters and initial conditions of Ackle and Ke "7 . Nonlinear reaction and nonlocal diffusion are given 

by 

a(^) :— max I £, — I + rno, for all ^, (5.5) 



/(w) — aw :— rw{K — w), for all w, (5.6) 

where s, mg, r and k are constant and positive parameters. We remark that in the numerical example 

of Ackle and Ke [2^ , the authors consider a nonlocal diffusion given by the expression — j . In our 

Jo u{x)dx 
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Figure 1: Extinction of the population density u{x,t) (left) and v{x,t) (right) at time T = 0.2. 



case the parameter e is a very small parameter and it plays a practical computational role to avoid the 
numerical overflow on the diffusion when the extinction of the population occurs, that is when w « or 
w ~ 0. We consider a parameter mo to the numerical study of an extinsion case of population and for a 
persistence case of population. In fact, the exponential decaying of the energy ()4.ip . can be interpreted 
as the extinsion of two populations u and v. That is occurs when the hypothesis m > Cp{Mi + a) of the 
Theorem 4.1, is verified. If m is too small, the decaying of the energy is not guarantee and a population 
persistence can be occur as we see in Figure 2. 

The initial condition is given by uo{x) — (5sin(7rx), with S = 1.95 (see [2]), and V(){x) = —uo{x) = 
—SsmijTx). In this example, a negative vq has not a physical or biological sense, if u and v represent 
densities of population, but we want to focus in the importance of the hypothesis of the Theorem 4.1, 
showing numerically that the exponential decaying of the energy does no occur when the hypothesis is 
not verified. We choose the parameter e ~ 10^^ for the nonlocal diffusion function (|5.5p . and we choose 
r = 1.0 and k — 10.0 for the reaction function (|5.6p . The discretization is given by J = lO'' and K — W^; 
we solve the linear system (|5.1l) - (|5.3p using Thomas algorithm programming in FortranOO. 

We consider here 2 simulations, the first one with mo = 1.0 (see Figure 1), and the second one with 
mo =0.1 (see Figure 2). The population persistence phcnomenom does not occur with a choice of a 
big amplitude 6 for the initial condition as it occurs in [2J. In fact we made simulations with different S 
values and the exponential decaying is always observed. 
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Figure 2: Persistence of the population density u{x,t) (left) and v{x,t) (right). 
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